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We study the two-point Topological Charge Density Correlator (TCDC) in lattice QCD with two 
degenerate flavours of naive Wilson fermions and unimproved Wilson gauge action at two values 
of lattice spacings and different volumes, for a range of quark masses. Configurations are gen- 
erated with DDHMC algorithm and smoothed with HYP smearing. In order to shed light on the 
mechanisms leading to the observed suppression of topological susceptibility with respect to the 
decreasing quark mass and decreasing volume, in this work, we carry out a detailed study of the 
two-point TCDC. We have shown that, (1) the TCDC is negative beyond a positive core and radius 
of the core shrinks as lattice spacing decreases, (2) as the volume decreases, the magnitude of the 
contact term and the radius of the positive core decrease and the magnitude of the negative peak 
increases resulting in the suppression of the topological susceptibility as the volume decreases, 
(3) the contact term and radius of the positive core decrease with decreasing quark mass at a 
given lattice spacing and the negative peak increases with decreasing quark mass resulting in the 
suppression of the topological susceptibility with decreasing quark mass, (4) increasing levels of 
smearing suppresses the contact term and the negative peak keeping the susceptibility intact and 
(5) both the contact term and the negative peak diverge in nonintegrable fashion as lattice spacing 
decreases. 
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1. Introduction 

As a consequence of the reflection positivity and the pseudoscalar nature of the relevant local 
operator in Euclidean quantum field theory, the two-point Topological Charge Density Correlator 
(TCDC) is negative at arbitrary non-zero distances [1]. In the continuum theory, close to the origin 
the two-point TCDC is negative and singular. From power counting, the singularity ~ — | _x; | s up to 
possible logarithms and hence is non-integrable. In order to obtain a positive and finite space-time 
integral (susceptibility), the TCDC should have a positive non-integrable singularity at the origin 
[1, 2]. However, it is possible to give a rigorous definition of topological susceptibility in Lattice 
QCD without power divergences using Ginsparg -Wilson fermion [3,4]. 

As the authors of Ref. [1] pointed out long time ago, this has non- trivial consequences for the 
derivation and interpretation of the Witten-Veneziano expression [5] for the 7]' mass. The negativity 
of the TCDC also has non-trivial consequences related to the nature of topological charge structure 
in QCD vacuum [6]. 

The issues related to two-point TCDC are best studied in the theory rigorously formulated on 
a Euclidean lattice. However, the lattice theory defined by a particular action may not be reflection 
positive. Fortunately, this is not a concern for the Wilson fermion. However, the breaking of 
chiral symmetry by Wilson term may lead to uncancelled divergences in topological susceptibility. 
Thus it is important to calculate topological susceptibility with Wilson fermion to check whether 
the cancellation indeed happens so that Wilson lattice QCD belongs to the same universality class 
as continuum QCD. The lattice operator for the topological charge density q(x) may extend over 
several lattice spacings, and thus for sufficiently small x, the continuum like behaviors are not 
expected. Nevertheless, continuum properties are expected to emerge as lattice spacings become 
smaller and smaller. Specifically, on a lattice with lattice spacing a, TCDC remains positive within 
a radius r c , which is expected to shrink to zero as a — > 0. The first investigation of lattice spacing 
dependence of the radius of the positive core and the negativity beyond the positive core of TCDC 
in lattice QCD was carried out in Ref. [7] in the context of overlap based topological charge density 
in quenched QCD. Later, similar study was carried out [8] for a variety of lattice QCD actions with 
and without quarks where discretization errors appear only at &{a ). 

Flavour singlet axial Ward-Takahashi identity relates the topological susceptibility %, which 
is the four-volume integral of TCDC, to the chiral condensate in the chiral limit [9, 10]. As a 
consequence, % vanishes linearly in the quark mass in the chiral limit. Furthermore, at a given 
value of the quark mass, % is suppressed as volume decreases [10, 11]. As part of an on-going 
program [12, 13] to study the chiral properties of Wilson lattice QCD (unimproved fermion and 
gauge actions), recently, we have demonstrated the suppression of topological susceptibility with 
decreasing quark mass in the case of unimproved Wilson fermion and gauge action [14, 15] where, 
the suppression of % with decreasing volume was also shown. In order to shed light on the mech- 
anisms leading to these suppressions, in this work, we carry out a detailed study of the two-point 
TCDC. 

2. Measurements 

We have generated ensembles of gauge configurations by means of DDHMC algorithm [16] 
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Table 1: Lattice parameters and simulation statistics. Here block, N2, N tr j and T refers to HMC block, 
step number for the force F 2 , number of HMC trajectories and the Molecular Dynamics trajectory length 
respectively. 

using unimproved Wilson fermion and gauge actions with rif ■ = 2 mass degenerate quark flavours. 
At p = 5.6 the lattice volumes are 16 3 x 32, 24 3 x 48 and 32 3 x 64 and the renormalized quark 
mass ranges between 25 to 125 MeV (MS scheme at 2 GeV). At j8 = 5.8 the lattice volume is 
32 3 x 64 and the renormalized physical quark mass ranges from 15 to 75 MeV. The lattice spacings 
are determined using nucleon mass to pion mass ratio and Sommer method. These determinations 
agree for the value of Sommer parameter ro = 0.44 fm. The lattice spacings at jS = 5.6 and 5.8 are 
0.069 and 0.053 fm respectively. The number of thermalized configurations ranges from 7000 to 
14000 and the number of measured configurations ranges from 200 to 500. 

The topological susceptibility % = Jd 4 x C(r) with the TCDC, C(r) = (q(x)q(0)), r = \x\ 
where q(x) is the topological charge density. For q(x), we use the lattice approximation developed 
for SU(2) [17], modified for St/ (3) [18] and implemented in the MILC code [19]. Unless otherwise 
stated we have used 3 smearing steps[20] in all our calculations. 

3. Results 

In order to extract the topological charge density reliably on the lattice, using the algebraic 
definition, smearing of link field is essential. Smearing however smoothens out short distance 
singularities. Excessive smearing may in fact wipe out the fine details of the singularity structure. 
Both the positive and negative contributions to % are affected in this manner. This is illustrated in 
Fig. 1 (left) where we show the effect of 3, 5, 10 and 20 HYP smearing steps on C(r) at jS = 5.8, 
K* = 0.15462 and lattice volume 32 3 x 64. However the susceptibility is remarkably stable under 
smearing after three smearing steps as illustrated in Fig. 1 (right) (taken from Ref. [15]). 

As was already stated, from theoretical considerations, we expect suppression of % with de- 
creasing volume at a fixed quark mass. In Fig. 2 (left) we present the finite volume dependence 
of the C(r) at j8 = 5.6 and K = 0.158 at lattice volumes 16 3 x 32, 24 3 x 48 and 32 3 x 64. We find 
that as volume decreases, the magnitude of the contact term and radius of the positive core de- 
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Figure 1: (left) Effect of smearing on C(r) at j3 = 5.8 and K = 0.15462 at lattice volume 32 3 x 64. (right) 
Effect of smearing on the % at j3 = 5.8 and K — 0.15462 at lattice volume 32 3 x 64 (taken from Ref. [15]). 
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Figure 2: (left) Finite volume dependence of the C(r) at j3 =5.6 and K = 0. 158 at lattice volumes 16 3 x 32, 
24 3 x 48 and 32 3 x 64. (right)Comparison of the radius of the positive core of C(r) at two different lattice 
spacings for comparable pion mass. Lattice volume is 32 3 x 64. 



crease and the magnitude of the negative peak increases resulting in the suppression of topological 
susceptibility as volume decreases. 

In Fig. 2 (right), we compare the radius of the positive core of C(r) at jS = 5.6 and 5.8 for 
comparable pion masses in physical units. The lattice volume is 32 3 x 64. The figure clearly 
exhibits the shrinking of the radius of the positive core of C(r) in physical units as one approaches 
the continuum. 

In order to understand the detailed mechanism behind the suppression of topological suscepti- 
bility with decreasing quark mass, we need to investigate the quark mass dependence of the various 
features of the C(r). In Fig. 3 (left) we present the quark mass dependence of C(r) with emphasis 
on the positive region at /3 = 5.6 and lattice volume 24 3 x 48. The magnitude of the contact term 
C(0) is seen to decrease with decreasing quark mass. In Fig. 3 (right) we present the quark mass 
dependence of C(r) with emphasis on the crossover from positive to the negative region of C(r) 
and the negative peak region at jS = 5.6 and lattice volume 24 3 x 48. The radius of the positive 
core and the magnitude of the negative peak of C(r) are seen to decrease and increase respectively 
with decreasing quark mass. The features presented in Figs. 3 result in the suppression of the 
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Figure 3: The quark mass dependence of C(r) with emphasis on the positive region (left) and with emphasis 
on the crossover from positive to the negative region of C(r) and the negative peak region (right) at j3 = 5.6 
and lattice volume 24 3 x 48. 
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Figure 4: (left) Topological susceptibility at j3 =5.6 and lattice volume 24 3 x 48 as a function of the 
quark mass mass, (right) Lattice spacing dependence of the contact term and the negative peak of C(r) at 
comparable pion mass for /3 = 5.6 and 5.8 and lattice volume 32 3 x 64. For comparison, the corresponding 
quantities for pure gauge lattice theory at j3 = 6.0983 (a = 0.078fm) and lattice volume 24 3 x 48 are also 
shown. 



topological susceptibility with decreasing quark mass. MILC collaboration [21] has made a simi- 
lar observation regarding the dependence of the negative peak on quark mass. In Fig. 4 (left) we 
present the corresponding topological susceptibilities (j8 = 5.6 and lattice volume 24 3 x 48) as a 
function of the quark mass which clearly shows the suppression as quark mass decreases. This 
figure includes K = 0. 158 in addition to the fc's presented in Figs. 3. 

In Fig. 4 (right) we show the lattice spacing dependence of the contact term and the negative 
peak of C(r) at comparable pion mass in physical units for j3 = 5.6 and 5.8 and lattice volume 
32 3 x 64. For comparison, the corresponding quantities for pure gauge lattice theory at /3 = 6.0983 
(a = 0.078fm) and lattice volume 24 3 x 48 are also shown. Both the contact term and the negative 
peak of C(r) increase with decreasing lattice spacing, in accordance with the expectation from the 
continuum theory. 

From the definition of the topological susceptibility % = J d A x C(r) = J 2n 2 {r i )dr C(r) it is 
useful to define [21] a local susceptibility 

X (r)= [ r 2x 2 (r' 3 )dr>C(r>) (3.1) 
Jo 
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Figure 5: The function %(r), defined in Eq. (3.1) as a function of r (left) and contributions to the suscepti- 
bility from positive and negative parts of C(r) (right) at j3 = 5.6 and 5.8 at comparable pion masses. 



in order to exhibit the lattice spacing dependence more clearly. In Fig. 5 (left) we plot %(r) versus 
r at two lattice spacings at comparable quark masses. The contribution from the positive part of 
C(r) results in a peak at short distance. This is followed by a decrease due to the negative part of 
C(r). As lattice spacing decreases, the contribution from the positive part increases resulting in the 
increase of the peak of x( r )- 

Define the contributions to the susceptibility from the positive and negative parts of C(r) as 
Xp = Jq c 2n 2 (r' 3 )dr J C(V) and Xn = X~ 2tt 2 (r /3 )^r r C(r') respectively. According to the expecta- 
tions from continuum theory, the negative singularity close to the origin and the positive singularity 
at the origin are both nonintegrable. Thus the contributions to x from positive and negative parts 
of C(r) are expected to diverge, nevertheless resulting in a finite x due to cancellation. In Fig. 5 
(right), we plot the contributions to the susceptibility from positive and negative parts of C{r) at 
/3 = 5.6 and 5.8 at comparable pion masses. The data exhibited in Fig. 5 (right) are in accordance 
with these expectations. 

It is known that the topological susceptibility decreases with decreasing quark mass and de- 
creasing volume. To understand the mechanisms leading to these suppressions, in this work, we 
carry out a detailed study of the two-point TCDC. We have shown that, with naive Wilson fermion 
and gauge action, (1) the two-point TCDC is negative beyond a positive core and radius of the core 
shrinks as lattice spacing decreases, (2) as volume decreases, the magnitude of the contact term and 
the radius of the positive core decrease and the magnitude of the negative peak increases resulting 
in the suppression of topological susceptibility as volume decreases, (3) the contact term and ra- 
dius of the positive core decrease with decreasing quark mass at a given lattice spacing and the 
negative peak increases with decreasing quark mass resulting in the suppression of the topological 
susceptibility with decreasing quark mass, (4) increasing levels of smearing suppresses the contact 
term and the negative peak keeping the susceptibility intact and (5) both the contact term and the 
negative peak diverge in nonintegrable fashion as lattice spacing decreases. 
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